Effective potential for classical field theories subject to stochastic noise 
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Classical field theories coupled to stochastic noise provide 
an extremely powerful tool for modeling phenomena as di- 
verse as turbulence, pattern-formation, and the structural de- 
velopment of the universe itself. In this Letter we sketch a 
general formalism that maps such systems into a field the- 
ory language, and demonstrate how to extract the one-loop 
physics for an arbitrary classical field theory coupled to Gaus- 
sian noise. The amplitude of the noise two-point function 
serves as the loop-counting parameter and is the analog of 
Planck's constant h. We define the effective action and the 
effective potential, and derive a general formula for the one- 
loop effective potential of a classical field theory coupled to 
translation- invariant additive Gaussian noise. 
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Classical field theories subject to stochastic noise [0,1), 
[e.g., stochastic partial differential equations (SPDEs)], 
provide powerful means for modeling a host of physi- 
cal phenomena ranging from turbulence |^], to pattern- 
formation [^, and to the structural development of the 
universe itself In this Letter we provide a brief 

summary of more extensive work to be presented in [^ . 
Consider the general class of SPDEs 



D(l){x, t) = F[(j){x,t)]+T^{x, t). 



(1) 



Here D is any linear differential operator, involving arbi- 
trary time and space derivatives, which does not explic- 
itly involve the field (j). The function F[(j)] is any forcing 
term, generally nonlinear in the field 0. The forcing term 
will typically not contain any time derivatives, but this is 



not an essential part of the analysis. The function 77 (af, t) 
denotes the source of noise. As yet the nature and proba- 
bility distribution of the noise need not be specified. We 
will take this SPDE and analyze it using functional in- 
tegral techniques. In this Letter we sketch out the field 
theory via the most direct route, avoiding the Martin- 
Siggia-Rose construction (with its extra unphysical con- 
jugate fields used for book-keeping purposes [00]). We 
will assume that the partial differential equation (|l|) , plus 
initial and boundary conditions, is well-posed. Given a 
particular realization of the noise, 77, the differential equa- 
tion is assumed to have a unique solution ipsoinix, t\r]). 
For any function Q{(j)) of the field (j) we can define the 
stochastic average, (over the noise), as 



(Q((/.)) EE / {V7J)V[V] Q(</.soln(f,t|??)). 



(2) 



Here V[ri] is the probability density functional of the 
noise. It is normalized to 1, but is otherwise completely 
arbitrary. We next use a functional delta-function to 
write the following identity 

0soln(^, = / (25'/') (j) S[(j) ~ (j)soln{x,t\ri)] 



[V^) 5[D^ - F[4>] - 77] ^JT^. (3) 
The Jacobian functional determinant, , is defined by 
J = det (d-^-^\. (4) 



The above is the functional analogue of a standard delta- 
function result. If f{x) = has a unique solution at a; = 
xq, then ccq = ^ dx x 5{x--xq) — J dx x S{f{x)) |/'(a;)| = 
J dx X S{f{x)) f'{x)[f'{x)]* . It is easy to see that one 
also has the identity 



Q{<t'soln{x,t\ri)) 



{V<jy) Q(0) 5[Dcj, - - ry] ^JJJ^ . 



The ensemble average over the noise becomes 



(5) 



(g(0)) = / {v^) VH x 

\v<i>) Qicj,) 5[D<f> - F[0] - 77] ^JJ^. (6) 
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The noise integral is easy to perform, with the result that 



(g(0)) = / iVcj,) V[Dcp - Q(0) x/jjt. (7) 



The effect of the noise appears in the stochastic average 
only through its probability distribution V[D(j} — F[(j)\]. 
The presence of the functional determinant is essential; it 
is often discarded without comment in the literature, but 
in general it must be kept to ensure proper counting of 
the solutions to the original SPDE. A particularly useful 
quantity is the characteristic functional. With an obvious 
condensed notation, dx = (V^x dt, we define 



Z[J] = ( exp / J{x) (j){x) dx 



J{V(t,)V[D(j)- F[^]]exp(^J J (j)dx^ VJ^- 



(8) 



When there is no risk of confusion we will suppress the dx 
completely. This key result will enable us to calculate the 
effective action and the effective potential in a direct way. 
At this stage we make some assumptions about the noise: 
We assume the noise is Gaussian with zero mean (so that 
the only non- vanishing cumulant is the second order one). 
If the noise has a non-zero mean, one can always redefine 
the forcing term F[(j)\ to make the noise be of mean zero. 
We still do not need to assume translation invariance, 
nor do we need to assume the noise is white, power-law, 
or colored. Arbitrary Gaussian noise is enough since it 
allows us to write the noise probability distribution as 



V[v] 



^det(27rG^) 
1 
2 



exp 



dx dy r]{x) G^{x,y) -qiy) 



(9) 



The characteristic function is thus 
1 



Z[J] 



yJdet{2TTGr, 
1 
2 



exp 



{Dq^ - F[cb])G-\Dq^ - F[cl,]) 



(10) 



This characteristic function contains all the physics de- 
scribable by (|l|). The noise rj has been completely elimi- 
nated and survives only through the explicit appearance 
of its two-point function Gri- Since the characteristic 
function is now given as a functional integral over the 
physical field 4>, all the standard machinery of statistical 
and quantum field theory can be brought to bear. This 
formula for the characteristic function demonstrates that 
(modulo Jacobian determinants) all of the physics of any 
stochastic differential equation can be extracted from a 



functional integral based on the "classical action" (note: 
this is not yet the effective action) 



{D<l>-F[4>])G-\Dc^-F[4>]) . (11) 



This action S[(j)\ can be viewed as a generalization of 
the Onsager-Machlup action ||l^ . Onsager and Machlup 
dealt with stochastic mechanics rather than field theory, 
and limited attention to white noise. As they developed 
their formalism with the notions of linear response theory 
in mind, they tacitly assumed the "forcing term" F[<j)\ to 
be linear, so that both the noise and the field fluctuations 
were Gaussian. 

The Jacobian functional determinant is sometimes 
(but not always) constant. This is a consequence of the 
causal structure of the theory as embodied in the fact 
that we are only interested in retarded Green functions. 
The situation here is in marked contrast to that in QFT, 
where the relativistic nature of the theory forces the use 
of Feynman Green functions {+ie prescription) In 
order to avoid enumerating special cases, and to have 
a formalism that can handle both constant and field- 
dependent Jacobian factors, we exponentiate the de- 
terminant via the introduction of Fadeev-Popov ghost 
fields 



1 



det (27r/) 



(V[g,g^]) exp 



D 



5F 



(12) 



The g and g'' fields are known as scalar ghost fields. They 
are defined in terms of anticommuting complex variables 
and behave in a manner similar to an ordinary complex 
scalar field except that there is an extra minus sign for 
each ghost loop. In all cases we are interested in the 
operator {D — SF/6(j)), even if not self-adjoint, is real and 
non-singular. This allows us to replace ^fj^ by J , so 
that the characteristic function (|l^) can be written as 



Z\J\ = 



1 
2 
1 
2 



i_/(P|0....-l)exp(/j0) 



exp 
exp 



{p<\>-F\c^\)G-\D<i>-F\<i}^) 
5F 



(13) 



This procedure trades off the functional determinants for 
two extra functional integrals. In order to develop a 
Feynman diagram expansion we treat the driving term 
F[<f>\ as the perturbation and expand around the free- 
field theory defined by setting F = 0. There are two 
propagators, one for the field, and one for the ghost 
fields. For simplicity, we will now take the noise to be 
translation invariant Translation- invariance lets us 
take simple Fourier transforms in the difference variable 
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X — y (more precisely: x — y and tx — ty), so that in 
momentum-frequency space we have 



D{-k,-uj) D{k,uj) 
1 



Gfioid(fc,t^) 

Gghost(fc,w) ^ 

D[k, u}) 

The Feynman diagram vertices are formally |^ 

D(fc,a;)0F[0] 



4> — F[(j3\ vertex: 
F[(j)\ - F[4>\ vertex: 
ghost vertex: 



Gr,{k,u) 
2 G^(fc, 



UJ 



1 ^sm 



'2^ 



(14) 
(15) 

(16) 
(17) 
(18) 



In practical calculations F[(j)] is a polynomial in the field 
and/or its gradients Some specific examples are dis- 
cussed in 13,1^]. 

To set up the effective action corresponding to (|^), 
and its loop expansion, it is useful to separate the two- 
point function for the noise into a shape, g2(x,y), and a 
constant amplitude, A 



, dcf 



G^{x,y) = Ag2ix,y). 



(19) 



We adopt the following simplifying normalization of the 
shape function g2{x,y) 



g^^{x,t) <:fxAt = 1 = g:^^ {k = {),uj = {)). (20) 



This is only a convention, not a restriction on the noise, 
since it only serves to give an absolute meaning to the 
normalization of the amplitude A Q . The virtue of in- 
troducing the amplitude parameter A is that it is the 
loop-counting parameter for the field theory discussed in 
this Letter To verify this, start by writing the charac- 
teristic function (with external sources rescaled for con- 
venience) 



1 



{'D[(t>,g,g'^]) exp 



^det[(27r)3G^] J " V ^ 

I j j {Dcj,~F[4>])g:,\D<l^-Fm 

1 
2 



exp 

exp 



g'iD-'-^^g 



(21) 



The generating function for connected correlation func- 
tions is defined by [lllJl^ 



W[J] ^ +A (In Z[J] - In Z[0]) , 



(22) 



and the effective action by |ll|l 



-W[J] 



J, ffliM^j, (23) 



with the stochastic expectation value of (j) for van- 
ishing external sources, J = 0. The previous equation 
defines the non-perturbative effective action. For spe- 
cific examples it is not very useful, and we often have to 
restrict ourselves to a perturbative calculation of the ef- 
fective action, after singling out an expansion parameter. 
One can always develop a Feynman diagram expansion 
provided that the classical action can be separated into 
a "quadratic piece" (to define propagators) and an "in- 
teracting term" (to define vertices). In this expansion 
the sum of all connected diagrams coupled to external 
sources J{x) is exactly I^[J], and the effective action 
r[(/);0o] corresponds to all (amputated) one-particle ir- 
reducible graphs (IPI), that is, Feynman diagrams that 
cannot be made disconnected by cutting only one prop- 
agator. 

To see that A is the loop-counting parameter, note 
that each field propagator is proportional to A while 
each ghost propagator is independent of A. The vertices 
that do not include ghosts are proportional to .4^^, while 
ghost vertices are independent of A. But each ghost loop 
contains exactly as many ghost propagators as ghost ver- 
tices. Thus if one assigns a factor A to each ghost prop- 
agator and a factor of A~^ to each ghost vertex, one will 
not change the total number of factors of A assigned to 
the Feynman diagram. Each diagram is proportional to 
A^~^ , where / is the number of (internal) propagators 
in the diagram and V is the number of vertices. 

It is a standard topological theorem that for any graph 
(not just any Feynman diagram) I—V = L—l, where L is 
the number of loops (l^Jl^ . This implies that stochastic 
field theories based on SPDEs exhibit exactly the same 
loop-counting properties as QFTs, except that the loop- 
counting parameter is now the amplitude of the noise two- 
point function (instead of Planck's constant h). The only 
subtle part of the argument has been in dealing with the 
Fadeev-Popov ghosts, and it is important to realize that 
this argument is completely independent of the details 
of the differential operator D and the forcing term F [</)]. 
When it comes to calculating the diagrams contribut- 
ing to the effective action, the extra explicit factor of A 



inserted in the definition of [equation (22)], guar- 

antees that the IPI graphs contribute to r[(/);(/)o] with a 
weight that is exactly A^ , which means that A is a bona 
fide expansion parameter. 

For any field theory, quantum or stochastic, the char- 
acteristic function can be written as 



Z[J] 



P0exp 



and the one-loop effective action as 

V[cj,; </)o] -5[0] + ialndet(52M) 
-(<^^</)o)-f 0(a2). 



(24) 



(25) 
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Here 52 5'^ S / 54>{x)54>{y) is the matrix of second func- 
tional derivatives of the action S[(j)\. The notation S[4>q\ 
is shorthand for S[{(t>[J = 0])]. When we consider field 
theories based on SPDEs, the loop-counting parameter a 
becomes A, and the bare action is replaced by M 



(26) 



The noise at this stage is translation-invariant and Gaus- 
sian. Inserting equation ( p6|) into the formula for the one- 
loop effective action, we obtain the general result (appli- 
cable to all SPDEs) 

+ 1^ {In det(52 [(/)]) - - In [</)]} 

-(0-^0o) + O(^2)^ (27) 

It is often useful to restrict to constant (spacetime inde- 
pendent) fields. For such field configurations the effective 
action reduces to the effective potential defined by 



, dof - 



(28) 



where il — J d'^x dt is the volume of spacetime. It is, in 
fact, sufficient to have both Dip = and — constant, 
the reason for this qualification will become clear when 
we discuss the effective potential for the KPZ equa- 
tion fill. A straightfor ward calculation ^ leads to 



= iF^[<f>] + ^A 



X In 



d'^fc duj 



i2n) 



d+l 



M + 0{A^). 



(29) 



The above formula contains the most important result of 
this Letter. This result is qualitatively very reminiscent 
of the effective potential for scalar QFT 



V[0] - V{cb) 



d'^k duj 



In 



(30) 



The major difference is the fact that the scalar QFT prop- 
agator is replaced with a more complicated and typically 
nonrelativistic propagator. Notice also that for SPDEs 
one can naturally adapt the noise to be both the source 
of fluctuations and the regulator to keep the Feynman 
diagram expansion finite. This follows immediately by 
inspection of ( p9| ) which shows that the (momentum and 
frequency dependent) noise shape function g2 will affect 
the momentum and frequency behaviour of the one-loop 
integral. The finiteness, divergence structure, and renor- 
malizability of this integral will very much depend on the 



functional form of §2 ■ Renormalizing such infinities on a 
case by case basis is relatively simple since we are only 
working to one-loop Jl3| , pt . 

In this Letter we have shown how to set up a "di- 
rect" functional formalism to explicitly calculate an ef- 
fective action and a corresponding effective potential for 
a general class of non-linear SPDE's with additive noise. 
The effective action and effective potential have the same 
structural definition as they do in quantum field the- 
ory, though here they can be calculated even when the 
SPDE does not follow from a Hamiltonian or Lagrangian 
variational principle. At the one-loop level in the noise 
amplitude we have derived a compact expression for the 
effective potential reminiscent of, but radically distinct 
from, the effective potential of scalar QFT. The effective 
potential for stochastic phenomena provides a powerful 
tool for exploring the influence of noise on a system, the 
onset of pattern formation, and noise-induced dynamical 
symmetry breaking. 

We hope to have convinced the reader that the "di- 
rect approach" described in this Letter, and developed 
more fully in is both useful and complementary to 
the more traditional MSR formalism Some ques- 

tions can more profitably be asked and answered in this 
"direct" formalism. For instance, our general formula 
for the one-loop effective potential for arbitrary SPDEs 
subject to translation-invariant Gaussian noise appears 
difficult to extract from the MSR formalism. 



Electronic mail: hochberg@laeff.esa.es 

Electronic mail: carmen@t6-serv.lanl.gov 

Electronic mail: mercader@laefr.esa.es 

Electronic mail: visser@kiwi.wustl.edu 
[1] P.C. Martin, E.D. Siggia, and H.A. Rose, Phys. Rev. 

A8 (1973) 423-508. 
[2] J. Zinn- Justin, Quantum field theory and critical pfie- 

nomena (Oxford University Press, Oxford, England, 

1996). 

[3] U. Frisch, Turbulence (Cambridge University Press, 

Cambridge, England, 1995). 
[4] M. Kardar, C. Parisi, and Y-C. Zhang, Phys. Rev. Lett. 

56 (1986) 889-892. 
[5] A. Berera and L-Z. Fang, Phys. Rev. Lett. 72 (1994) 

458-461. 

[6] D. Hochberg and J. Perez-Mercader, Cen. Rel. Gravit. 
28 (1996) 1427. 

[7] T. Coldman, D. Hochberg, R. Laflamme, and J. Perez- 
Mercader, Phys. Lett. A222 (1996) 177. 

[8] J.F. Barbero, A. Dommguez, T. Coldman, and J. 
Perez-Mercader, Europhys. Lett. 38 (1997), 637-642. 

[9] D. Hochberg, C. Molina-Pan's, J. Perez-Mercader, and 
M. Visser, Effective action for stoctiastic partial dif- 
ferential equations, Phys. Rev. E60 (1999) 6343-6360 



4 



| ;ond-mat/99042iq ]. 
[10] L. Onsager and S. Machlup, Phys. Rev. 91, 1505 
(1953). 

[11] S. Weinberg, The quantum theory of fields /,// (Cam- 
bridge University Press, Cambridge, England, 1996). 

[12] R. J. Rivers, Path integral methods in quantum field 
theory (Cambridge University Press, Cambridge, Eng- 
land, 1987). 

[13] D. Hochberg, C. Molina-Pan's, J. Perez-Mercader, 
and M. Visser, Effective potential for the massless 



KPZ system Physica A280 (2000) 437-455 |cond- 
mat/990439l|. 



[14] D. Hochberg, C. Molina-Pan's, J. Perez-Mercader, and 
M. Visser, Effective potential for the reaction- diffusion- 



decay system J. Statist. Phys. 99 (2000) 903-941 [cond- 
mat/9909078| . 



5 



